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ABSTRACT 

Near-infrared (2.2 jim) observations of a sample of 48 elliptical galaxies in the Coma 
cluster have been carried out and used to study the near-ir fundamental plane (FP) 
' of ellipticals in this cluster. An rms scatter of 0.072 dex was found for this relation, 

similar to that of its optical counterpart, using the same sample of galaxies. This 
corresponds to an uncertainty of 18% in distances to individual galaxies derived from 
' this relation. The sensitivity of the near-infrared FP to the star formation or changes 

\ in metallicity and stellar population among the ellipticals was explored and found to 

■ be small. Although, a likely source of scatter in this relation is contributions from the 

^i: ' Asymptotic Giant Branch (AGB) population to their near-infrared light. Allowing for 

, observational uncertainties, we find an intrinsic scatter of 0.060 dex in the near-ir FP. 

' The cluster galaxies presented here, provide the zero-point for the peculiar velocity 

I studies, using the near-infrared FP. 

. Changes in the slopes of the D — a and L — a relations of ellipticals between 

the optical and near-infrared wavelengths were investigated and found to be due to 
variations in metallicity or age (or a combination of them). However, it was not possible 
to disentangle the effects of age and metallicity in these relations. 
P ■ We find M/L oc M" with a = 0.18 ±0.01 at the near-ir and a = 0.23 ±0.01 at the 

optical wavelengths, using the same sample of galaxies. This relation is interpreted 
c/3 ■ as due to a mass-metallicity effect or changes in age or the IMF slope with mass. 



Using evolutionary population synthesis models, we find that the effects of age and 
metallicity decouple on the {M/L)k vs. Mg2 and {M/L)k vs. (V— K) diagrams. The 
models suggest that the observed trends on these relations may be due to an age 
sequence while metallicity mainly contributing to the scatter. 
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1 INTRODUCTION 

Elliptical galaxies define a Fundamental Plane (FP) on the 
3-parameter space of the half-light radius Re, the mean sur- 
face brightness within that radius {SB)e, and the central 
velocity dispersion a (Djorgovski & Davis 1987; Dressier et 
al. 1987a). The scatter around the FP is only ~ 0.08 dex in 
logi?e for cluster ellipticals. The existence of this scaling rela- 
tion provides one of the most important constraints to model 
the formation and evolution of elliptical galaxies (Bender et 
al 1993; Guzman et al. 1993; Renzini & Ciotti 1993; Zepf & 
Silk 1996). Also, it can be used as a distance indicator to 
trace the non-Hubble motions in the local universe (Dressier 
et al 1987b; Lucey & Carter 1988; Lynden-Bell et al. 1988), 



and potentially as a cosmological evolutionary probe when 
studying the FP for cluster ellipticals at higher redshifts 
(Franx et al. 1996; van Dokkum & Franx 1996; Barger et al. 
1996). 

The biparametric nature of elliptical galaxies suggests 
that the virial theorem is the main constraint on their struc- 
ture (Faber et al. 1987; Djorgovski 1987; but see Guzman et 
al. 1993). The tilt of the FP relative to the virial equation 
implies a systematic variation either in the IMF, the shape of 
the light profile, the age and metallicity or the dark and lu- 
minous matter distribution with galaxy mass (Guzman et al. 
1993; Ciotti et al. 1996). Our understanding of the physical 
mechanisms responsible for such variation is still at a very 
early stage. Moreover, substantial fine-tuning is required to 
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produce the tilt while preserving the small scatter of the 
galaxy distribution along the FP (Ciotti et al. 1996). The 
tightness of the FP for cluster ellipticals is clear evidence 
for a very standardized and synchronized production of el- 
lipticals, with the vast majority of their stellar populations 
being formed at redshifts a > 2, at least for those ellipti- 
cals in clusters (Renzini 1995; Ellis et al. 1996). However, 
there is growing evidence that ellipticals outside the core 
of rich clusters tend to have a younger stellar component 
(Rose 1985; Pickles 1985; Bower et al. 1990). As a result, 
non-cluster ellipticals at a given Re and/or a have, on aver- 
age, bluer colours, higher effective surface brightnesses, and 
lower Mg2 line strength indices as compared to their cluster 
counterparts (Larson et al. 1980; Guzman et al. 1992; de 
Carvalho & Djorgovski 1992). This systematic variation of 
stellar population properties due to the environment may in 
turn translate into a zero-point offset of the FP-based dis- 
tance indicators (such as _D„ — a) and thus lead to spurious 
peculiar velocities (Kaiser 1988; Silk 1989; dc Carvalho & 
Djorgovski 1992; but see Burstein et al. 1990). Indeed, using 
stellar population evolutionary models, Guzman & Lucey 
(1993) have combined the FP with the a — Mg^ relation 
to create a new distance indicator that is independent of 
age/environmental differences. When the new distance in- 
dicator is applied to the ellipticals with a younger stellar 
component, the large positive peculiar motions previously 
reported arc greatly removed (Guzman & Lucey 1993). 

Most studies on the FP have been done at optical wave- 
lengths. Observations in the infrared may shed light on the 
origin of the tilt and scatter of the FP by probing a wave- 
length region with a different sensitivity to stellar popula- 
tion effects. In particular, the effect of age and/or metallic- 
ity variations in ellipticals is expected to be greatly reduced 
in the K-band (i.e., 2.2 fim) since the galaxy light at these 
wavelengths is mainly sensitive to the old stellar compo- 
nent and is less affected by line-blanketing. Moreover, the 
K-band is significantly less sensitive to both the presence of 
dust in ellipticals and galactic extinction. The new K-band 
FP can also be used as an independent distance indicator 
to check for stellar-population effects on the peculiar veloc- 
ity measurements derived using the FP in the optical. Early 
work on the infrared FP relied on single element photom- 
etry and hence, lacked detailed surface photometry (Gregg 
1995a; Recillas-Cruz et al. 1990, 1991). Recently, Pahre et 
al. (1996) have shown that the coefRcients and scatter of the 
FP in the K-band for a sample of ellipticals in five nearby 
clusters are similar to those derived in the optical. They 
conclude that the observed small differences are consistent 
with the reduction of metallicity effects in the near-infared 
bandpass. 

In this paper, we present the first results of a large-scale 
K-band survey of elliptical galaxies located in different en- 
vironments (i.e., clusters, field and Great Attractor region). 
Here we study the coefficients and scatter of the K-band 
FP, using accurate surface photometry for a homogeneous 
sample of 48 galaxies in the Coma cluster. By considering 
a single cluster, we avoid introducing artificial scatter due 
to distance errors or non-universality of the FP. The aim of 
this study is three fold: 1) to shed light on the origin of the 
FP by using a wavelength region with a different sensitivity 
to the stellar population effects than the optical; 2) to es- 
tablish the intrinsic scatter on the K-band FP; and 3) to set 



the zero-point of the new infrared distance scale using the 
K-band FP for Coma cluster ellipticals. 

In the next section, the observations and data reduction 
are discussed. This is followed by the study of the near- 
infrared scaling laws in section 3. The implications for the 
origin of the FP and the new distance scale are discussed in 
section 4. Finally, our conclusions are summarized in section 
5. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Observations 

The observations used in this study were made at the United 

Kingdom Infrared Telescope (UKIRT) in the period 28-30th 
April 1994. Near-infrared K-band (2.2 /im) images of 48 el- 
liptical galaxies in the Coma cluster were obtained, using 
the 256 X 256 infrared array detector (IRCAM3) with a pixel 
size of 0.286". The galaxies are selected to be ellipticals (i.e. 
no lenticulars) from the sample in Lucey et al (1991, Table 
6) with reliable velocity dispersion and optical photometry. 
They are located at difi'erent distances from the center of the 
Coma cluster (a = 12 57 19 5 = 28 15 51 (1950)). Also, they 
have optical half-light diameters smaller than 1', fitting the 
field of view of the IRCAM3. An exposure time of 10 sec. 
was used with 6 co-adds. For each object, nine exposures 
were taken at different positions, separated by 20", giving a 
total integration time of 9 mins. per object. This technique 
was employed to avoid the effect of bad pixels and cosmic 
rays, to construct reliable fiatfields for each frame and to 
measure the sky background. Surface photometry of smaller 
galaxies can be affected by the seeing condition which was 
estimated to be around 1". The sensitivity of the results 
to this effect will be explored in the following analysis. Ob- 
servations of faint UKIRT standards were also carried out 
during the night and were used to monitor the accuracy of 
the photometry. 



2.2 Data Reductions. 

Observations were dark subtracted using the dark frames 
closest (in time) to each object. For each galaxy, the nine 
mosaic frames were then median filtered and the result nor- 
malised to its median pixel value to construct the flatfield. 
The dark subtracted, flatfielded frames were then registered 
and mosaiced to increase the S/N. The data reduction was 
carried out, using the CCDPACK software in the STAR- 
LINK environment. The median sky value was then esti- 
mated for each frame individually and used to sky-subtract 
that frame. The mosaic frames, providing a larger sky cov- 
erage than single frames, were used for this purpose and the 
sky values were measured as far away from the center of 
galaxy as possible. The variations in sky estimates at differ- 
ent positions on the frame were monitored, giving a measure 
of the uncertainties in these values. This did not exceed 0.02 
mag., even for the bigger galax;ies for which the uncertain- 
ties are larger. Therefore, this is taken as our estimated error 
due to sky subtraction. 

The standard stars were reduced in the same way and 
used to establish the extinction relation for each night. 
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Changes in the zero-points, estimated from the standards 
throughout the run (internal photometric error), is 0.02 mag. 

2.3 Photometry 

Aperture photometry was performed on the final images and 
used to construct the near-ir growth curves for individual 

galaxies. The magnitudes wore corrected for atmospheric ex- 
tinction. Galactic absorption (taking Ak = 0.085^b (Palirc 
et al. 1996) and the median Ab for the Coma galaxies from 
Faber et al (1989)) and redshift effect (assuming k{z) = 3.3z 
for ellipticals at the K-band (Persson et al (1979)). The pho- 
tometric accuracy is checked by comparing magnitudes of 
the objects in our sample with the same galaxies from other 
independent observations. The K-band magnitudes for 18 
galaxies in the present sample, in common with Bower et 
al (1992), are measured over an aperture of 17" diameter, 
used by these authors. The two magnitudes are compared 
in Figure la and show a mean offset of 

{Kthia study — KBower et al.) = —0.06 ± 0.01 

with an rms scatter of 0.03 mag. Assuming equal errors in 
both the data sets, this corresponds to an internal error 
of 0.03/v'2 = 0.02 mag. per gala^xy, consistent with that 
claimed by Bower et al (1992). The cause of the zero-point 
offset between the two data sets is not clear but is likely to 
be due to different detectors used. However, since the data 
here are 2-D and Bower et al's is based on single element 
photometry, we expect to have a more accurate control over 
the sky subtraction in our data. Neverthless, the offset is 
relatively small. 

The growth curves, corrected for Galactic absorption 
and redshift effects, were used to construct near-infrared 
surface brightness profiles for each galaxy, corrected for the 
(l-|-z)~* surface brightness dimming. Near-infrared isopho- 
tal diameters, corresponding to a mean surface brightness 
of 16.5 mag/arcsec"^ are then estimated. This isophote was 
adopted to be directly comparable with the mean optical 
surface brightness of 19.8 mag/arcsec.^ (assuming V — K = 
3.3 mag. for ellipticals) used in the optical Dy — cr rela- 
tion. To examine the accuracy of our surface photometry, 
we have also obtained data for a sub-set of 12 ellipticals in 
our sample, using the IRIS at the Anglo-Australian Tele- 
scope (March 1995). These diameters, measured at a mean 
K-band surface brightness of 16.5 mag./axcsec^, are com- 
pared with the present measurements in Figure lb and show 
a mean offset of 

{/o,9(D^^^^^/D^'^^)) = 0.013 dea; ± 0.006dea- 

with an rms scatter of 0.019 dex. Assuming that the uncer- 
tainty in these measurements is shared equally between both 
the data sets, we estimate an observational error of 0.01 dex 
in log(Z)if ) values. 

The diameters of smaller galaxies in our sample might 
be affected by the seeing condition. To estimate this effect, 
we used the models in Lucey et al (1991) and estimated the 
effect of 1 arcsec seeing on the aperture photometry (the 
observations here were done in 0.8-1.2 arcsec seeing condi- 
tions) . The correction was then applied to one of our smaller 
galaxies and its photometric parameters were recalculated. 
We find that at 1 arcsec seeing conditions, an increase in 
log(i3if) of only 0.01 dex was needed due to seeing. The 



diameters of larger galaxies were almost insensitive to this 
effect. The range in surface brightness values between 16 and 
17 mag/arcsec^ was explored in 0.1 intervals. The results in 
the following sections were found to be insensitive to the 
actual choice of the isophote. 

The near-infrared effective diameter and surface bright- 
ness (i.e. the diameter containing half the total luminosity 
of a galaxy and the mean surface brightness within that di- 
ameter) were calculated for galaxies in the present sample 
by fitting the observed K-band curve of growths to the dc 
Vaucouleurs r^^* law. The accuracy of this fitting proce- 
dure was estimated by generating simulated profiles follow- 
ing the r'^''* law with known effective parameters and artifi- 
cially added noise, analogous to the photometric errors. An 
r^^^ fit to these profiles recovers both the input effective pa- 
rameters, giving errors of 0.02 dex and 0.03 mag/arcsec^ in 
the effective diameter and surface brightness values respec- 
tively. The total K-band magnitudes are then estimated by 
extrapolating the r^''''-law fit for individual galaxies to large 
radii, using their corresponding values of the effective diam- 
eter and surface brightness. These magnitudes are corrected 
for Galactic absorption and redshift, using the relations dis- 
cussed above. Also, the effective surface brightness estimates 
are corrected for the (1 -I- a)"'* dimming effect. 

Finally, the K-band magnitudes, measured over an 
aperture of 20 arcsec. diameter are estimated and combined 
with their corresponding V-band data from literature, mea- 
sured over the same aperture (Lucey et al 1991). The V — K 
colours were then estimated and corrected for Galactic ex- 
tinction and redshift, using the above prescription for the 
near-ir and the relations from Faber et al (1989) for the 
optical data. 

2.4 The Catalogue 

The sample of 48 elliptical galaxies in the Coma cluster, 
observed in this study, is presented in Table 1. Column 2 
gives the distance (in degrees) of galaxies from the center of 
the cluster. In columns 3 and 4, the near-infrared (K-band) 
isophotal diameters in arcsec. (corresponding to a mean sur- 
face brightness of {SBk) ~ 16.5 mag/arcsec^) and total 
magnitudes of ellipticals are presented. These are corrected 
for redshift effect, surface brightness dimming and Galactic 
absorption. Column 5 gives the optical V-band diameters 
(in arcsec), corresponding to a mean surface brightness of 
(SBv) = 19.8 mag/arcsec^, taken from Lucey et al. (1991). 
The velocity dispersion of ellipticals (log(f7)- in Km/sec), 
listed in column 6, are taken from Lucey et al. (1997). For 
galaxies with no available velocity dispersion measurements 
from this source, the values in Lucey et al. (1991) were con- 
verted to the above system using the common galaxies be- 
tween the two samples. The velocity dispersions have an 
Eissociated error of 0.03 dex. Columns 7 and 8 present the 
near-infrared effective diameters (in arcsec.) and effective 
mean surface brightnesses. This is followed in columns 9 and 
10 by the K-band aperture magnitudes and V— K colours re- 
spectively, measured over am aperture of 20 arcsec. diameter 
and corrected for Galactic extinction and redshift. The Mg2 
line strengths are listed in column 11. These are also taken 
from Lucey et al. (1997). For galaxies with no such measure- 
ments from this reference. The Mg2 indices from Lucey et 
al. (1991) are used after conversion to this system. 
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Table 1. The near-infrared catalogue for elliptical galaxies in Coma 
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N4789 


1.525 


1.341 


9.30 


1.365 


2.416 


1.349 


16.53 


10.20 


3.25 


0.304 


N4971 


1.656 


1.080 


10.57 


1.072 


2.250 


1.096 


16.54 


11.21 


3.28 


0.291 


D159-89 


1.986 


1.023 


10.70 


1.075 


2.232 


1.159 


16.98 


11.37 


3.10 


0.273 


D159-83 


2.500 


1.182 


10.23 


1.152 


2.306 


1.088 


16.16 


10.82 


3.38 


0.275 


D160-159 


2.790 


1.062 


10.66 


1.073 


2.358 


1.089 


16.60 


11.22 


3.24 


0.280 


N4673 


3.294 


1.329 


9.68 


1.356 


2.345 


1.069 


15.51 


10.25 


3.18 


0.270 


D159-43 


4.545 


1.122 


10.60 


1.118 


2.399 


0.968 


15.93 


11.07 


3.33 


0.338 


D159-41 


4.750 


1.012 


11.11 


1.025 


2.277 


0.876 


15.98 


11.51 


3.23 


0.324 



3 THE NEAR-INFRARED SCALING LAWS 

3.1 The Near-Infrcired Fundamental Plane 

The near-infrared fundamental plane of ellipticals in the 
Coma cluster is established by fitting a plane to the dis- 
tribution of galaxies on the effective diameter (Ae), effective 
mean surface brightness {{SBK)e) and velocity dispersion 
(cr) space, using the data from Table 1. A simultaneous 3- 
parameter fit, using 48 galaxies, gives 

log{Ae) = (1.38 ± 0.26)log{a) + (0.30 ± Om){SBK}e + ci 

with an rms scatter of 0.072 dex in log(Ae). An edge-on 
view of the near-ir FP is shown in Figure 2. The shape of 



the near-infrared FP and its scatter is in good agreement 
with that found by Pahre et al (1996), using 12 ellipticals 
in the Coma cluster. The presence of curvature on the in- 
frared FP is explored by investigating the dependence of the 
scatter in Figure 2 on the FP parameters. No relation has 
been found between the residuals around the FP (A (FP) 
= 1.38 log{(T) + 0.3Q{SBK)e + ci-log{Ae)) and the effective 
diameter or surface brightness. The observed rms scatter in 
the noar-infrarcd FP corresponds to an uncertainty of 18% 
in distances to individual galaxies from this relation. 

To investigate differences between the infrared and op- 
tical FPs, we have also constructed the optical (V-bamd) 
FP, using the same sample of 48 ellipticals listed in Table 1. 
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The coefficients of the optical FP, estimated by performing 
a 3-parameter fit to these data, correspond to 1.44± 0.04 
and 0.32± 0.01 for log(a) and {SBv}e respectively. The rms 
scatter in the optical FP here is 0.074 dex, which is simi- 
lar to the rms scatter of 0.08 dcx (Djorgovski and Davics 
1987; Jorgensen et al 1996) and 0.07 dex (de Carvalho and 
Djorgovski 1992) for the B-baiid and 0.075 dex (Lucey et al 
1991) for the V-band FP, using other independent samples. 

Differences between the FP coefficients from different 
studies, found in literature, are mainly due to the adopted 
fitting mathods (i.e. taking Ae or a as the independent vari- 
ables in the least-squares solutions) or sample selection, with 
consideration to minimise the bias in the fit. For example, it 
has been estimated that a least-squares fit in log{Ae), pro- 
duces a bias of the order of 5% for the log{a) coefficient, 
caused by observational errors (Jorgenson et al 1996). In 
this study, it is attempted to reduce such biases by perform- 
ing a 3 parameter fit to the FP. However, adopting different 
fitting methods, the FP coefficients hero lie in the range 
1.32-1.50 {log{a)) and 0.30-0.32 (< SBe >) at the near- 
infrared and 1.23-1.44 ilog{a)) and 0.30-0.35 (< SB^ >) 
at the optical wavelengths. This shows the sensitivity of the 
FP coefficients to the adopted fitting technique, with the 
values in the above range, being consistent with the present 
study. The FP parameters found here, also agree closely with 
other independent estimates of the near-infrared (Palire et 
al 1996) and optical (Faber et al 1987; Guzman et al 1993) 
FPs. 

The coefficients of the near-ir and V-band FPs are sim- 
ilar despite different sensitivities to the stellar population 
and metallicity at these wavelengths. However, a smaller 
rms. scatter might have been expected in the near-ir FP 
since at this wavelength, the light samples a more uniform 
stellar population and is less affected by differences in line 
blanketing among the ellipticals. These will be addressed 
later in this section. 

Studies of the optical FP have revealed the presence 
of a young stellar population in ellipticals to be partly re- 
sponsible for the observed scatter in this relation (Gregg 
1995b). Moreover, at a given a, Coma ellipticals with lower 
Mg2 values have, on average, slightly larger Dv diameters, 
with most of them located outside the core of the cluster 
(i.e. r>l deg.), indicating contributions from an intermedi- 
ate age stellar population, preferentially in ellipticals in low 
density environments (Guzman et al 1992). This confirms 
that Mg2 features arc sensitive to both the 'young' stel- 
lar population and local environment of ellipticals with the 
galaxies having larger contributions from the 'young' stel- 
lar populations satisfying 5.2Mg2 - log(Dv) — 0.442 < -0.2 
(Guzman and Lucey 1993). There arc only five galaxies in 
our near-infrared sample which satisfy this relation (N4876, 
N4807,N4789, 159-83 and N4673) with four of them located 
at the outskirts of the cluster (r> 1 deg.). These galaxies 
are indistinguishable from the rest of the ellipticals on the 
near-ir FP in the Coma cluster (Fig 2). 

The sensitivity of the near-ir FP to changes in the stel- 
lar population and metallicity among the ellipticals are fur- 
ther investigated by comparing the A (FP) values with the 
Mg2 index and V-K colour residuals (at a given a) from 
the Mg2 — a and (V-K)— cr relations respectively (Figure 3). 
No relation has been found. A similar study of the residu- 
als diagram, using the optical FP consisting of both field 



and cluster ellipticals, shows a significant trend which is 
mainly interpreted as a stellar population efi'ect (Prugniel 
and Simien 1996). However, using the sample of cluster ellip- 
ticals in this study, we find no relation between the residuals 
from the V-band FP and A (Mg2) or A (V— K) estimates. 
This result has been confirmed independently by Jorgensen 
et al. (1996), using a different sample of elliptical galaxies in 
clusters. Therefore, it is not clear if the lack of correlation on 
the residuals diagrams in Figures 3a and 3b, using the near- 
infrared FP in the Coma cluster, is due to the uniformity of 
our cluster galaxies or, represents a genuine effect, implying 
that the observed scatter on this relation is not affected by 
metallicity or contributions from the young stellar popula- 
tion. These scenarios arc currently explored by extending 
the study of the near-infrared FP to the field ellipticals. 

However, it is likely that changes in contributions from 
the intermediate age Asymptotic Giant Branch (AGB) stars 
to the near-infrared light of ellipticals is partly responsible 
for the observed scatter in this relation. Spectroscopic ob- 
servations of ncar-ir CO (2.3 fim) absorption features of el- 
lipticals in clusters have revealed these objects to have a 
smaller population of AGB stars than those in general field 
(Mobashcr and .Tames 1996). Therefore, it may be possible 
to introduce ncar-ir CO features as an extra parameter to re- 
duce the scatter on the infrared FP in the same way the Mg2 
indices were introduced to the optical FP to correct the relar- 
tion for contributions from the young population (Guzman 
and Lucey 1993). 

Considering the measurement errors of 0.02 dex in 
log(Ae), 0.03 dex in log((T), 0.03 mag/arcsec^. in the {SBK)e 
and 0.02 mag. due to sky subtraction, and considering the 
fact that the errors in log(Ae) and {SBK)e are not indepen- 
dent (Kormcndy 1977), we estimate an observational error 
of 0.040 dex for the near-ir FP. Compared with a total rms 
error of 0.072 dex, this gives an intrinsic scatter of 0.060 dex 
in the near-ir FP. 



3.2 The Dk - CT Relation 

The Dk — cr relation for our sample of 48 ellipticals in the 
Coma cluster is shown in Figure 4. Taking a as the inde- 
pendent variable (Lyndon-Bell et al. 1988), a least squares 
fit to the Dk — cr relation, using all the 48 galaxies, gives 

log{DK) = (1.37 ± 0.08)log((T) - (2.09 ± 0.16) 

with an rms scatter of 0.073 dex. To explore the wavelength 

dependence of this relation, the optical Dy — o" relation, 
using the same sample of 48 ellipticals from Table 1, is con- 
structed. This is consistent with a slope of 1.15 ± 0.09 and 
rms scatter of 0.077 dex, in agreement with studies based 
on other independent samples (Jorgensen et al 1996; Lucey 
et al 1991; Djorgovski and Davics 1987). 

The similarity in the rms scatter between the two rela- 
tions implies that dynamical parameters or the presence of 
an intermediate-age stellar population contributing to the 
near-ir light of ellipticals are responsible for the observed 
scatter in the Dk ~ f relation. 

The slight wavelength dependence of the slope of the 
isophotal diameter-velocity dispersion relation (~ 2a effect) 
leads to a correlation between log{Dv/DK) and log{a)- 
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(Figure 5) as 

log{Dv/DK) = (-0.18 ± Om)log{a) + (0.42 ± 0.05) 

This implies that over the range of velocity disper- 
sions covered in our sample {A(loga) ^ 0.6), the diame- 
ters change by 10% (ie. more massive galaxies have corre- 
spondingly larger near- infrared diameters). Using the stel- 
lar synthesis models of Worthey (1994), the slope of the 
log(Z)y/Dx)— log((T) relation, due to changes in metallicity 
and age, is predicted and compared with the observed re- 
lation in Fig 5. The models correspond to a range in age 
from 5 to 17 Gyr at a constant metallicity of [Fe/H]=0.25 
and variations in metallicity from —0.25 to 0.5 at a fixed 
age of 12 Gyr. The two models arc normalised to have 
log(Dy/£>x) = -0.01 at 12 Gyr and [Fe/H]=0.25. These 
normalisations are chosen to be close to the values expected 
for local ellipticals. The predicted relations are consistent 
with the observed slope, implying that changes in age (ie. 
stellar population), metallicity or a combination of them 
among the ellipticals is responsible for the wavelength de- 
pendence of the slope of their diameter-velocity dispersion 
relation. However, it is not possible to disentangle the effects 
of age and metallicity using these models. 

3.3 The Lx - cr Relation 

The total near-infrared luminosities of elliptical galaxies in 
the Coma cluster are used from Table 1 to establish the 
near-ir Lk — cr relation. The total magnitudes axe employed 
here in order to have this relation on the same magnitude 
scale as the previous studies of its optical counterpart. The 
Lk — o relation is presented in Figure 6. A least squares fit, 
taking a as the independent variable and using 48 ellipticals 
in the Coma, gives 

Ktot = (-6.78 ± 0.50)log(a) + (26.28 ± 1.14) 

with an rms scatter of 0.50 mag. Using the aperture mag- 
nitudes (measured over am aperture of 20 arcsec. diameter- 
column 4 of Table 1) gives a shallower slope of —5.35 ± 0.37 
and a smaller rms scatter of 0.32 mag. The optical rela- 
tion, using the total V magnitudes for the same sample of 
ellipticals has a slope and rms scatter of —6.51 ± 0.57 and 
0.49 mag. respectively. The similarity in the slope and rms 
scatter between the near-infrared and optical L — o rela- 
tions implies that dynamical and structural characteristics 
in ellipticals are probably as important as the photometric 
parameters (i.e. stellar population, metallicity and age) in 
defining the mass-luminosity relation in these systems. This 
will be explored in the following discussion. 

Taking the observational errors in K (0.03 mag.), 
Zo5(cr)-(0.03 dex) and sky background estimate (0.02 mag.), 
gives an observational error of 0.19 mag. in the Lk — o re- 
lation. Compared with the observed scatter of 0.50 mag., 
an intrinsic scatter of 0.46 mag. in the Lk — c relation is 
estimated. 

It is known that the colours of elliptical galaxies are 
correlated with their absolute magnitudes and indirectly, 
with their mass, in both optical (Faber 1973; Visvanathan 
1983) and infrared (Persson ot al 1979; Mobasher et al 1986) 
wavelengths. The existence of such a relation has impor- 
tant implications for studies of evolution of these galax- 
ies. This colour-absolute magnitude (mass) correlation can 



therefore be translated into a more direct relation between 
the colours and velocity dispersions in ellipticals. Using the 
V— K colours measured over an aperture of 20 arcsec. diam- 
eter for the sample of 47 galaxies with available such data 
from Table 1, the {y — K) — log{a) relation is presented in 
Figure 7. A least squares fit gives 

V-K = (0.74 ± 0.10)log{a) + (1.51 ± 0.19) 

with an rms scatter of 0.083 mag. in V—K colours. Using the 
stellar synthesis models of Worthey (1994), the {V — K) — a 
relations due to changes in age and metallicity are predicted 
and compared with the observed relation in Figure 7. The 
models are calculated for a given metallicity, [Fo/H]=0.25, 
and a change in age from 5 to 17 Gyrs, and a range in 
metallicity,— 0.25 < [Fe/H] < 0.5, at a constant age of 12 
Gyrs. The models are normalised to have [Fe/H]=0.25 and 
an age of 12 Gyrs at V — K = 3.35 mag., the values expected 
for local ellipticals. They are consistent with the observed 
(V — K) — cr relation, again indicating that changes in either 
age or metallicity (or a combination of them) among the 
ellipticals are, at least, partly responsible for this relation. 
However, these models cannot disentangle the effects of age 
and metallicity on the {V — K) — a plane (Fig. 7). Recent 
studies also propose the presence of a mass-metallicity rela- 
tion among the cluster ellipticals (Kodama & Arimoto 1996) 
with the metallicity being the more fundamental parameter. 

To explore if the above results are model dependent, 
we investigate the behaviour of the observable parameters 
in ellipticals to changes in age and metallicity, using in- 
dependent stellar synthesis models (Bruzual and Chariot 
1996). At a fixed age of 12 Gyrs and over the metallic- 
ity range —1.64 < [Fe/H] < 1.008, we predict a slope of 
A{y - K)/A{Mg2) = 5.23 on the (V-K)-Mg2 plane which 
is close to A{V — K)/A{Mg2) = 6 found for a fixed metal- 
licity of [Fe/H]=0.093 and a range in age from 5 to 12 Gyrs. 
These results are used to predict the A{V — K) / A{log{a)) 
slope (ie. the relation in Fig. 7), using the empirical Mg2- 
<7 relation derived from the present sample. We find that 
the (V— K)-(T relations, predicted by Bruzual and Chariot 
(1996) model, have similar slopes due to changes in metal- 
licity (at a given age) or age (at a given metallicity). This 
confirms that the results in Figures 5 and 7 are independent 
of the models used (ie. the metallicity and age trends are 
degenerate) . 

The {V — K) — a relation is independent of distance and 
hence, can be used to study evolution of ellipticals at dif- 
ferent environments (i.e. in field and clusters) and redshifts. 
This relation can also be used to estimate velocity disper- 
sion of ellipticals at high redshifts (for which a spectroscopic 
measurement of a is difficult), using their V— K colours and 
allowing for evolution with redshift. The accuracy in log(a') 
values, estimated from the (V— K)-cr relation, is 0.11 dex. 

3.4 The Near-infrared M/L Ratios 

The near-infrared FP found here differs, at a high confi- 
dence level, from that expected from a pure virial theorem 
if light directly traces the mass and if ellipticals arc dy- 
namically and structurally homologous (f e oc a^I^^; where 
Ie is the effective surface brightness in units of luminos- 
ity per area and is the effective diameter). This leads 



© 0000 RAS, MNRAS 000, 000-000 



The Near- Infrared Fundamental Plane of Ellipticals 7 



to a relation between the near-ir mass-to-luminosity ratios 
log{M/L) = 2loga - log{D^/2) + QA{SB)e + oi and mass 
log{M) = log{De/2) + 2log{a)+a2 of the ellipticals (Fig. 8)- 
(Faber et al 1987; Djorgovski 1987). The {M/L)k and mass 
estimates in Figure 8 arc in Solar units, assuming z = 0.023 
for the Coma and Ho = 50 km/sec/Mpc. This relation is 
consistent with {M/L)k oc M" where a = 0.18 ± 0.01, in 
excellent agreement with a = 0.16 ± 0.01 estimated from a 
K-band study of 59 ellipticals in 5 neaxby clusters (Palire et 
al 1996; see also Recillas-Cruz et al 1990, 1991). 

To investigate the origin of this relation and explore 
if it is caused by differences in the dark-to-luminous mat- 
ter distribution among the ellipticals or is due to changes in 
metallicity, age or stellar population, we extend this study to 
optical wavelengths, using the same sample of galaxies (Fig. 
8). The optical relation has a steeper slope of a = 0.23±0.01, 
indicating that over the same range in mass, the optical and 
infrared M/L ratios change by factors of 3 and 2.3 respec- 
tively. The slight wavelength dependence of the slope of the 
M/L vs. M relation implies that changes in the luminous- 
to-dark matter distribution and/or luminosity profile among 
the ellipticals is not the only parameter responsible for the 
observed relations. Other observable parameters (see below), 
are also likely to be as important. 

In order to have an entirely independent measurement 
of the slope of the M/L vs. M relations in the two pass- 
bands, we estimate both the M/L ratios and masses of in- 
dividual galaxies, using their respective effective diameters 
in the ncar-ir and optical wavelengths. However, this intro- 
duces the unphysical effect of the dependence of the dy- 
namical mass on the wavelength in Figure 8. Therefore, the 
shift in the masses of the same object between the optical 
and near- infrared (M/L) vs. M relations in Figure 8 is due 
to differences in the shape of the surface brightness profiles 
in the two bands, and to observational uncertainties in the 
values of the effective diameters. Nevertheless, using the ef- 
fective diameters measured only in one passband to estimate 
the mass of galaxies, the slope of the M/L vs. M relations 
remain the same as above, indicating that the conclusions 
in this section do not depend on different effective diameters 
for individual galaxies, used to estimate the mass. 

The effects of stellar population, metallicity and age on 
the near-infrared mass-to-luminosity ratios of the ellipticals 
are investigated in Fig 9, using V-K colours and Mg2 line 
indices. Employing the stellar synthesis models of Worthcy 
(1994), changes in the {M/L)k with age and metallicity are 
calculated for 5 Gyrs < t < 17 Gyrs at constant metallicities 
of [Fe/H]=-0.25, and 0.25 and for a range in metallicity 
corresponding to —0.5 < [Fe/H] < 0.5 at a constant age 
of 12 Gyrs. The models, normalised to have log(M / L) k = 
—0.5 at 12 Gyrs and [Fe/H]=0, are compared with the data 
in Figs 9a and 9b. 

There is only a weak dependence of the noar-ir mass- 
to-luminosity ratio on metallicity. Furthermore, the effects 
of age and metallicity in ellipticals, as shown in Fig. 9, 
decouple when using their (M/L)k estimates (ie. the age 
and metallicity sequences are orthogonal). Assuming the 
same 'collapse factor' (the ratio of the radius covering half 
the total mass to the effective radius) for all the ellipti- 
cal galaxies in our sample. Figures 9a and 9b suggest that 
the age is responsible for the observed trends on both the 



(M/L)K-Mg2 and (M/L)if-(V-K) diagrams while metal- 
licity mainly contributes to the scatter in these relations. 



4 IMPLICATIONS OF THE NEAR-INFRARED 
FUNDAMENTAL PLANE 

The near-infrared FP is less affected by differences in line 
blanketing, metallicity and stellar population among the el- 
lipticals, compared to its optical counterpart. This implies a 
FP relation with a potentially reduced scatter at the near- 
infrared wavelength. Therefore, the similarity of the rms 
scatter between the near-ir (0.072 dcx) and optical (0.074 
dex) FPs found here, was unexpected and could have im- 
portant implications for constraining models of formation of 
ellipticals. Understanding the origin of this scatter provides 
a challenge in studying the evolution of elliptical galaxies. 

Study of the scatter around the optical FP with elliptic- 
ity (ie. shape) has found no relation between the two quan- 
tities (Jorgensen et al 1996; Jorgensen et al 1993). As dis- 
cussed in section 3.4, it is likely that changes in the {M/L)k 
among ellipticals is responsible for at least part of the scatter 
in their ucar-infrarcd FP. This may be caused by differences 
in age and/or metallicity (Figure 9), contributions from the 
Asymptotic Giant Branch (AGB) population to their near- 
infrared light or differences in matter distribution and the 
internal dynamics (ie. orbital anisotropy or rotation; Bender 
et al 1993) among the ellipticals. However, similarity of the 
scatter between the optical and near-infrared FP relations 
indicates that structural and dynamical differences rather 
than changes in agc/motallicity or the AGB contributions 
among the ellipticals arc probably responsible for the ob- 
served scatter in their near-infrared FP. 

The tightness of the FP has been exploited in estimat- 
ing relative distances between elliptical galaxies in the field 
or clusters. However, the optical FP is affected by contribu- 
tions from the young population or residual star formation. 
Indeed, it is likely that the relatively blue ellipticals, found 
in the Great Attractor region, axe responsible for the pro- 
posed streaming motion (Guzman and Lucey 1993). The 
near-infrared FP may be loss sensitive to this effect (see sec- 
tion 3.1) and hence, the FP at this wavelength is expected 
to have a more stable slope and zero-point. Therefore, this 
relation can be used to test the reality of the streaming 
motion and its dependence on the non-uniformity of stel- 
lar populations among the ellipticals. For this purpose, we 
have performed near-ir (K-band) imaging observations of a 
sub-sample of ellipticals in the Great Attractor region. Us- 
ing the present sample of ellipticals in the Coma cluster, we 
establish the zero-point of the field near-ir FP in order to 
use this to explore the sensitivity of the streaming motion to 
changes in stellar populations among the ellipticals. These 
results will be presented in a forthcoming paper. 

Recent studies have indicated that elliptical galaxies 
are not the homogeneous population once assumed and are 
likely to have a complex evolutionary history involving in- 
teraction, starburst and infall of material (White and Frenk 
1991; Cole et al. 1994). This leads to changes in the lumi- 
nosity, mass and morphology of these galaxies with time, 
implying that their progenitors may be distinctly different 
from the galaxies we see today. A study of the FP of ellip- 
ticals at high redshift is needed to address these scenarios. 
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Recently, using the high resolution imaging capability 
of the Hubble Space Telescope and spectroscopic facilities 
on the largest ground-based telescopes (MMT and Keck), 
it has become possible to explore the evolution of the op- 
tical FP and M/L of ellipticals with redshift (Franx et al 
1996; van Dokkum & Franx 1996; Barger et al 1996). This 
reveals a FP similar to those at low redshift with a compa- 
rable scatter. However, the observed change in the optical 
M/L is found to be smaller than that expected from pop- 
ulation synthesis models based on a given metallicity, IMF 
and formation redshift (van Dokkum & Franx 1996). With 
the installation of the near-infrared (NICMOS) detectors on 
the HST, it is now possible to extend this study to longer 
wavelengths. The near-infrared light mainly measures the 
integrated star formation and hence, is more closely related 
to the mass function of the ellipticals. Therefore, study of 
the evolution of near- infrared FP and M/L of ellipticals with 
redshift constrains the infall and merging in these galaxies 
and is needed for a careful interpretation of the evolution of 
the luminosity function of galaxies. 



5 CONCLUSIONS 

The results of this study can be summarised as follows: 

(i) the near-infrared (2.2 ^.m) fundamental plane of ellip- 
tical galaxies in the Coma cluster is studied, using a sample 
of 48 galaxies. This shows an rms scatter of 0.072 dex, simi- 
lar to its optical counterpart, corresponding to an accuracy 
of 18% in distances derived from this relation. 

No relation is found between the scatter around the near- 
infrared FP and the V— K colours or Mg2 line strengths in 
ellipticals, probably implying negligible contribution to this 
relation due to changes in metallicity or the young popula- 
tion among the Coma cluster ellipticals. The AGB stars are 
proposed as the likely source responsible for the observed 
scatter in the near-ir FP. Allowing for observational uncer- 
tainties, we estimate an intrinsic scatter of 0.060 dex in this 
relation. 

(ii) The Dk — u and Lk — u relations are studied. These 
have an rms scatter of 0.073 dex and 0.50 mag. respectively, 
compared to scatters of 0.077 dex and 0.49 mag. for their op- 
tical counterparts. The wavelength dependence of the slopes 
of these relations lead to correlations between {Dv / Dk) and 
(V-K) with the mass (cr) of ellipticals. Using stellar synthesis 
models, we interprete these as due to an age or metallicity 
(or a combination of these two) effect. However, it is not pos- 
sible to disentangle the effect of age and metallicity, using 
these parameters. 

(iii) The relation between the M/L ratio and the mass 
of ellipticals [M/L oc M") in the Coma cluster is studied. 
The slight difference in the slope of this relation between 
the infrared {a = 0.18 ± 0.01) and optical (q = 0.23 ± 0.01) 
wavelengths is caused by a mass-metallicity relation and not 
due to changes in the luminous-to-dark matter distribution 
among the ellipticals. 

(iv) Using the stellar synthesis models, we disentangle 
the effects of age and metallicity on the (M/L)k— Mg2 and 
(M/L)x — (V— K) diagrams. The trend in these relations is 
interpreted as an age sequence while metallicity mainly con- 
tributs to its scatter. We propose the use of (M/L)_r- to sep- 
arate the effects of age and metallicity in elliptical galaxies. 



(v) We propose to use the near-infrared FP of ellipticals 
in the Coma cluster to fix the zero-point of its corresponding 
field sample in the Great Attractor region. This is required 
in order to explore the effect of stellar population and metal- 
licity on the estimates of the streaming motion of galaxies. 
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6 FIGURE CAPTIONS 

Figure 1. (a) Comparison between the K-band magnitudes 
of the Coma ellipticals in common with galaxies observed in 
Bower et al (1992). The magnitudes are measured over an 
aperture of 17" diameter, used in Bower et al. 
Figure 1. (b) Comparison between the K-band isophotal 
diameters (at 16.5 mag/arcsec.^) of a sub-set of 12 ellipticals 
in the Coma cluster observed both at the UKIRT (IRCAM3) 
and the A AT (IRIS). 

Figure 2). An edge-on view of the ncar-infrarod fundamen- 
tal plane for 48 ellipticals in the Coma cluster. The line is 
a plane fit to the data. The rms scatter in log(Ae) is 0.072 
dex. 

Figure 3. (a). Deviations from the near-infrared funda- 
mental plane {A{FP) = 1.38 log{a) + 0.30 < SBk >e 
+ci — log{Ae)) are plotted against A (V— K) colour residuals 
from the (V— K)-log(cr) relation. The lack of correlation im- 
plies that metallicity or age do not significantly contribute 
to the obseved scatter on the near-infrared FP of ellipticals 
in the Coma cluster. 

Figure 3. (b). The same as Fig. 3a for the A (Mg2) resid- 
uals from the Mg2-log((T) relation. Again, lack of correla- 
tion indicates that changes in the young population or pos- 
sible environmental dependence among the ellipticals in the 
Coma cluster do not contribute to the scatter on their near- 
infrared FP. 

Figure 4). The log(Dif )-log(cT) relation for 48 ellipticals in 
the Coma cluster. The T)k diameters are measured at 16.5 

mag/arcsec.'^ isophote. 

Figure 5). The log(Dv/DK )-log(cT) relation for ellipticals 
in the Coma cluster. The lines are calculated from the stel- 
lar synthesis models of worthey (1994). They correspond to 
changes in age (from 5 to 12 Gyrs) at a constant metallic- 
ity of [Fe/H]=0.25 (dotted line) and changes in metallicity 
(-0.25 < [Fe/H]< 0.5) at a constant age of 12 Gyrs (sohd 
line). The models are transformed to the observed parame- 



ters, using the empirical log((7)-Mg2 and log(Dy/DK)- (V- 
K) relations (from table 1), with the Mg2 and V-K values for 
a given age and metallicity taken from Worthey (1994) mod- 
els. The models are normalised to have log(D\//Dx)=— 0.01 
at [Fc/H]=0.25 and 12 Gyr age. 

Figure 6). The Ktot— log(cr) relation for the Coma cluster 
ellipticals. The total K magnitudes are used. The relation 

has an rms error of 0.5 mag. 

Figure 7). The (V— K)— log(a) relation for the Coma ellip- 
ticals. Colours arc measured over an aperture of 20" diame- 
ter. The lines are calculated from the stellar synthesis mod- 
els of Worthey (1994). They correspond to changes in age 
from 5 to 17 Gyrs at a constant metallicity of [Fc/H]=0.25 
(dotted line) and changes in metallicity from [Fe/H]=-0.25 
to [Fe/H]=0.5 at a constant age of 12 Gyr (solid line). The 
model estimates for a are calculated using the empirical re- 
lation as in Fig. 5. It is not possible to disentangle the effects 
of age and metallicity from this relation. 
Figure 8). Correlation between the near-infrared (•) and 
optical (*) mass-to-luminosity ratio; log (M/L) — 2loga — 
log {Dj2)+0A(SBe}+ai and mass; log (M) = log{De/2) + 
2log{a) + 02 for the Coma ellipticals. The same sample of 
galaxies are used in both cases. The zero-points correspond 
to ai = —10.67 and —11.24 for the near-infrared and optical 
wavelengths respectively and a2 — 5.47. Both the mass-to- 
luminosity ratio and mass are in solar units. For each wave- 
length, the mass corresponding to the effective diameter for 
that wavelength is plotted, leading to a slight difference in 
the mass estimated for the same objects in the near-infrared 
and optical relations. 

Figure 9 (a). Changes in the near- infrared mass-to- 
luminosity ratios in Solar units (calculated as in Fig. 8) with 
the V— K colours. The lines arc predictions from the stel- 
lar synthesis model of Worthey (1994) and correspond to a 
change in age from 5 to 17 Gyrs at a constant metallicity of 
[Fe/H]=-0.25 (dashed hnc) and [Fc/H]=0 (dotted line) and 
a change in metallicity from [Fc/H]=— 0.5 to [Fe/H]=0.25 
at a constant age of 12 Gyrs (solid line). The models are 
normalised to have log(M/L)if =— 0.50 at [Fe/H]=0 and an 
age of 12 Gyrs. 

Figure 9 (b). The same as Fig 9a for the Mg2 features. 
The lines are predictions from the stellar synthesis model of 
Worthey (1994) and correspond to a change in age from 5 
to 17 Gyrs at a constant metallicity of [Fe/H]=0.25 (dashed 
line) and [Fe/H]=0 (dotted line) and a change in metal- 
licity from [Fc/H]=— 0.5 to [Fe/H]= 0.5 at a constant age 
of 12 Gyrs (solid line). The models are normalised to have 
log(M/L)K=-0.50 at [Fe/H]=0 and an age of 12 Gyrs. 
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